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Abstract 

We study long time dynamical properties of a chain of harmonically bound 
Brownian particles. This chain is allowed to wander everywhere in the plane. 
We show that the scaling variables for the occupation times Tj, areas Aj and 
winding angles 9j (j = 1, . . . , n labels the particles) take the same general form 
as in the usual Brownian motion. We also compute the asymptotic joint laws 
P({Tj}), P({Aj}), P({0j}) and discuss the correlations occuring in those distri- 
butions. 



1 Introduction 

In order to study the dynamics of dilute polymer solutions, P.E. Rouse proposed in 1953 
his famous model of harmonically bound Brownian particles (Rouse chain) [1]. Since 
that time, this model has become very popular in the field of polymer science. It appears 
that, despite its drawbacks and limitations (in particular, absence of excluded volume 
and hydrodynamic interactions), it is conceptually important and useful to study the 
dynamics of polymers in melts [2,3]. In this paper, we will consider the free planar 
motion of such a chain of n particles (monomers) and especially address its long time 
(t — > oo) properties from the Brownian motion viewpoint. 

A configuration of this chain being represented by a complex n-vector z (the compo- 
nents Zi, i — 1, . . . , n are the complex coordinates of the particles), we will study closed 
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trajectories of length t, i.e. z(t) = z(0) or open ones (z(0) fixed, z(t) left unspecified, 
i.e. integrated over). 

More precisely, if we consider some given bounded domain S of area S and define 
the occupation time Tj as the time spent inside S by the j th particle, our goal is to 
compute the joint probability distribution P(Ti,T 2 , . . . ,T n ) (= P({7}})). Similarly, A,- 
and 9j being respectively the area enclosed by the trajectory of the j th particle and its 
winding angle around O, we will be interested in the joint laws P({Aj}) and P({9j}). 

On general grounds, we expect that the various properties of the chain will be strongly 
influenced by the free Brownian motion of the center of mass (c.o.m.) of the chain. 
But will they strictly satisfy the same laws ? With the same scaling variables ? And 
what about the correlations among the different variables ? Before answering those 
questions, we first recall some standard results concerning a planar Brownian particle 
with a diffusion constant D [4-8] . 

Results i) and ii) concern open trajectories when t — > oo while iii) concerns closed 
trajectories and is valid for all times t: 

i) Kallianpur-Robbins' law [4] for the probability distribution of the occupation time 
T of a bounded domain of area S: 

ii) Spitzer's law [5] for the angle 9 wound around a given point: 
with the characteristic function: 

(e^') = e-l A l (3) 

iii) Levy's law [6] for the area A enclosed by the closed trajectory of the particle: 

P { A ~ 2Dt) ~ 2 coshVA') ^ 

iBA'\ _ (t) /-v 

The distributions i) and iii) have moments of all orders in contrast with ii) that has 
none. 

Those laws were discovered more than 40 years ago and since that time, many refine- 
ments have been made. For instance, in [7], the authors found the asymptotic (t — > oo) 



joint law of the small and big (#+) windings. 0_( resp. 0+) are the angles wound 
around and only counted when r is smaller (resp. greater) than some fixed r (r is 
the distance separating the particle from O). With the rescaled angles 6' ± = ^y, the 
characteristic function writes [7]: 

ii { x-e>_ + x + e> + )\ = 1 (6) 

X 1 cosh(A+) + ^sinh(A+) 

(A+ = A_ = A gives back Spitzer's law (2)). 

Remark that (2) and (6) don't depend on the diffusion constant. This is quite 
different from the Brownian motion on a bounded domain surrounding O. In that case, 
we have [9]: 



>^ xe ) = e - cD l A l* (7) 

where c is a constant depending on the geometry and the boundary conditions. Here, 
D, appears as a multiplicator of |A|. We will use this remark at the end of the paper. 



2 The free Rouse chain 



Starting our study, we consider the following set of coupled Langevin equations: 



z\ = k (z 2 - zi) + 771 

z x = k (z l+1 + - 2z t ) + rji , 2 < I < n - 1 (8) 

Zn = k (z n -i — Z n ) + T) n 

where k is the spring constant and r\ m (= i] mx + ir) my ) a gaussian white noise: 



(Vm(t)) = 

(r] m (t)Vm'(t')) = 28 mmf 5{t-t') (9) 

(This noise would correspond to a D — 1/2 diffusion constant if particles were free). 

For the chain c.o.m., we get z G = \ (T,7=iVi) = Vg with (r) G (t)r) G (t')) = Hit - t'). 
The c.o.m. motion is free with D = l/(2n). 

Introducing the complex n- vector 77, eq.(8) can be written in a matrix form: 



z — — k~M z + T] 
where M is the tridiagonal (n x n) matrix: 
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with eigenvalues: 



2 1 - cos 



<3 ~ 1) 
n 



(wi = 0; det' M = n™ =2 = n ) 

With the matrix u = diag(ajj), we can write: 

uj = R 1 M R 

z = R Z 



(11) 



(12) 
(13) 

where R is an orthogonal matrix and the components of Z are the normal coordinates, 
that we will widely use in the sequel. From Rji = j = 1, . . . ,n, we deduce that 
Zi (= E"=i Zi/y/n) is essentially the c.o.m. coordinate. Remark also that E™=2 w i|^| 2 — 

Er= 2 ki-^-ii 2 = 'zMz. 

Let us call V(z, z( \t) the probability for the chain to go from configuration z^ at 
t = to z at time t. P satisfies a Fokker-Planck equation [10]: 



d t V = [% k M z +* d g k M z + 2 *^ <9 2 J P (14) 

where <9 Z (resp. d z ) is a n- vector of components d Zi (resp. <9 f J and t d z (resp. t d z ) is 
the transpose of d z (resp. <9 2 ). The solution can be written in terms of a path integral 
(VzVz = X$ =1 VziDzi): 



V(z,z^,t) 



with 



det (e* fcM ) £ o VzVzexp (-^ £ \z + kMz)(z + kMz) drj (15) 
i^z< >).G (*,* (0) ,f) 



G (z,z(°\t) 



e -k(tg Mz _t z (o) Mz (<0) = 

e -iEr= 2 (i^-ii a -i^ 0) -^ii 2 ) = e -!Er =2 ^(i^i 2 -i^ (0) i 2 ) 

jT } £>z£>,2 exp (-^ £ ('I i + k 2 *z M 2 z - 2k TrM) dr) = 



= z 



-tff 



,(0) 



-2 { ^9 z + -r z M 2; — /c TrM 
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(16) 
(17) 



In fact, V, eq.(15), can be easily deduced from the gaussian distribution of rj (use 
(10); det(e tfcM ) is the functional Jacobian for the change of variable r\ — > z [11]). 

Go(z, z^ \ t) is most conveniently written in terms of the normal coordinates Zi and 
zf\ clearly exhibiting the free motion of the c.o.m. [12]: 



G (z,z^\t) = ^e-^~ z ^ 2 x 
Zirt 

x ft exp [-\ (z tCl Z t + zf Ci zf - zf Si Z { - Z iSi zf) }) (18) 

; Ci = kuJi coih.(kuL>it) 



sinh(kuJit) 

When kt 3> 1, we get, for G , the limiting expression: 



G?(z,z<»,t) = ^ e ^i^-^i 2 nf^ e -¥(i^i 2 + i^i 2 )') (19) 

= g (z,z^,t).g (z,z^) 
where Q is the c.o.m. propagator and g can be simply written in terms of the zf 



g (z,z {U >) =n\-\ e V J 



(20) 



Furthermore, as can be easily checked, V is properly normalized: 
JdzdzV(z,z(°\t) = 1 (G given by (18) or (19)). 

Now, we turn to the computation of the joint law P({Tj}). 

3 Occupation times distribution 

Recall that Tj is the time spent by particle j inside a bounded domain S of area S. We 
consider trajectories starting at t — from some given configuration z^ and reaching 
at time t the final configuration z. Leaving z unspecified, we have, with positive pi 's: 



e-£r =1 ^)=det (e tkM ) x 



X 



J dzdz j* VzVzexp (- J \z + kMz)(z + kMz) + V P (z)^j drj (21) 
= J dzdz F(z, z {0) ) G P (z, z^\ t) (22) 
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with 



G P (z,z^,t) = (z 



-t(H +V P ) 



AO) 



Vp(z)=I>M*) 



(23) 
(24) 



ls(zi) is the indicatrix function of the domain S. Symbolically, we write: 



Gp= E(- 1 ) mG o(VpG )' 

m=0 



(25) 



with 



G (V P G ) 



ft fU, 

/ dt m / 

Jo Jo 



X 



dt m _i • • • y o " d*i / ^11 d^d^J G (*, ^ (m) , t - t ro ) 
x ^(zHjG^H^f-D, tm _ ^^(-D) . . . ^ (1) )G (* ( V (0) , *i) (26) 

is the chain configuration at time t,-; dz^dz^ = Yl^dz^dz^). 

Let us compute the contribution N m (t) of this generic term to (22). Integrating over 
z, we have: 



Aj) 



X 



ft rt m rt 2 r I m \ 

N m (t) = (~i) m j o dt m j o dt m ^...j o dt,] n* 1 ^ w ^ (o) ) 

x V P (z^)G (z^\ z^~ l \ t m - t m _!) . . . V P (zW)Go(zW, z®\ h) (27) 
= (-l) m £dt m J ^fi^^^Vi'^)^ .PC^, «W)*(t m , {^*>}) (28) 

$ is a time convolution product of the free propagators Gq. Disregarding for the moment 
the spatial integrations, the above expression is well-suited, in the limit t — > oo, for 
applying Tauberian theorems [13]. Introducing the Laplace Transform $: 

$(u,{z®})= e- ut <l>{t',{z^})dt' =Y[G {z^,z^ 1 \u) (29) 
• /() A i 

we notice that, when u — > + : 



Go(^ (fc) ,^ (fc_1) ,«)= / e-^Go(^ fc) ,2 (fc - 1) ,Od* / ~ / e-^Go(^ (fc) ,^ (fc_1) ,Od* / (30) 

JO 
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for some large a. So, we can use the asymptotic form in the computation of Go and 
get: 

Go(s<*\ z^\u) ~_ + In (I) i- 5o (^ (fc) , (31) 
A weak Tauberian theorem [13] gives for the time integration in (28): 

jf dt m $(t m ,{z«}) Jlg (zW,z< k -V) (32) 

Finally, the result for N m (t) is: 

AUO (-1)- (=^)" (33) 

iF = /fnd^- i j Vf-wfj") (34) 

Lp is computed with Vp, eq.(24): Lp = (I]" =1 p\)<S 
Rescaling the occupation times 7]: 

r < = ££i (35) 

we get: 

oo / n \ m i 

This relationship is actually valid, by analytic continuation, for all the positive pi 's (and 
not only when J2?i < !)■ This is because the distribution P({X,}) has moments of all 
orders and consequently (e~^=i PiTi ) is holomorphic when Re(pj) > 0. 

(36) leads to the probability distribution: 

n 

P({T1}) = 6{T[) J] Sm - TU) (37) 

i=2 

(n — 1 gives back the Kallianpur-Robbins' law). 

So, in the large time limit, the T,- L 's are strongly correlated leading to identical (T/) 
's. Moreover, we remark that Tj scales like n and, also, that the law for the c.o.m. would 
be the same as for one monomer (compare (35) to - (1) with D = l/(2n)): the c.o.m. 
free motion dominates this process. 

We also got similar exponential distributions for the rescaled variables T 1 in the 
following cases: 
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i) T is the time spent when the whole chain is inside S. Lp, eq.(34), is now computed 
with Vp(z) = p (niLi ls( z i))- Introducing 

W ( S ) = j (jldzidzMz^ e -*Er=2l**-^-il a (38) 

the rescaled variable writes: 

r ry 1 2 * T (39) 

" \k) nw(S) hit 1 j 

In contrast with (35), k is now present in the asymptotic law. For instance, if S is 
a small disk of radius r (kr^ -C 1), then w(S) ~ S n and T scales like k 11 ^ 1 : when 
k grows, the chain collapses and it is easier to confine it inside a given domain. 

ii) T is the time spent by the chain when at least one of its particles is inside S. T' 
is similar to (39) except that w(S) must be changed: Lp is now computed with 
V P (z)^p(l-U2 =1 (l-ls(zi))). ' 



To conclude this section, let us draw two lessons: 

i) The scaling variables take the same general form as for the free Brownian particle. 
In the sequel, we will show that it is still true for the other quantities we study. 

ii) For the computation of the perturbation theory, when t — > oo, we can system- 
atically use the asymptotic form of the unperturbed propagator. Obviously, 
for this consideration to hold, we must be sure that the perturbation series is well 
behaved. In those conditions, we will make a wide use of this remark. 



4 Areas distribution 

Now, we compute the areas distribution P({A,}) for closed trajectories of length t 
starting and ending at some fixed z^°\ To do so, we insert the constraint: 




in the measure (15) and use the relationship 5(x) — ^ / e tBx dB. It is easy to show that 
this manipulation amounts to add n different magnetic fields Bj to the initial system. 
Those fields are uniform, orthogonal to the motion plane and such that particle j is 
submitted to Bj. 
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With the (n x n) diagonal matrix B (B^ = Bi5ij), we get 



, ; i 2tt j \G (zW,zW,t), 
with G B (z<-°\z<-°\t) = ( | e~ tHB | z<® ) (42) 
# B = #o + Vb (43) 

V B (*) = Jf-^B^+^-B^) + i^B 2 ^ (44) 
2 ^ ' 8 

= -(- t ZB'd z + t ZB'd 2 ) + - t ZB" Z (45) 

(B' = R ^BR, B" = R^R). In principle, P({Aj}) depends on z^ but we will show 
that, actually, this is not the case when t — > oo (remark that, for all t, P({Aj}) doesn't 
depend on the c.o.m. of z^: this is due to translation invariance and this is the reason 
why we consider the propagator and not the partition function that would diverge like 
the area of the plane, leading to serious problems in the perturbation theory). 

Now, let us sketch the perturbative computation of Gb- Following our previous 
remarks, we will use G™ for the unperturbed propagator. The generic term writes: 

(-l) m J*dt m J* m dt rn _ 1 ...J* 2 dt 1 J (j[dl&)dz® 
. . . G™(z( j+1 \ z®, t j+1 - t 3 )VB{z^)G™{z^\ z^ x \t 3 - . . . (46) 

With normal coordinates, Vb(z^) takes the form: 

M^) = £ Q {-Z?B' u d zlj) + zfB' a d^) + \zTKZ?) (47) 

We now proceed by inspection: 

i) For the terms Z^B' u d v (j), only will survive, after integration, those contributions 

(j\ 

with % = j. The same holds for Z, B'^d-u). Moreover, the diagonal contributions 



exactly cancel except when % — j — 1. 



ii) We reach the same conclusion for the terms Z^B" iX z\^ (non diagonal terms vanish 
after integration. Diagonal contributions, i — j > 1, are subleading - compared 
to i — j — 1 - when kt >> 1: we recover the fact that the process is dominated by 
the c.o.m. motion). 

We are finally left with an effective perturbation V-^: 

Vb S = \ {-ZiB' u d Zl + ZxB'n^J + ^B^lZxl 2 (48) 
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Only the first mode is affected by the magnetic fields and we can disregard the other 
modes that will cancel when taking the ratio Gb/G . 

Remark that 

B ;i =I(g B( ) m dB ?1 = I(g B ?) (49) 

The effective hamiltonian for the remaining mode writes: 

H* = -2^% i + i (-ZxB'nfcx + ^iB'n^ + ^(B' n ) 2 ^i| 2 ) +\ (B" n - (B' n ) 2 ) |Z X | 2 

(50) 

It describes the behavior of a charged particle sub mitted to an u niform magnetic field 
B' n and an harmonic oscillator of frequency uj = \\J'Q'{ 1 — (B' n ) 2 . Using known results 
about this problem [14], we immediately get: 



G B (* (0) ,* (( V) ty/B» u 
G (zW,zW,t) ' ^ ' — " 




y/B^ (cosh - cosh (§B' n )) \Z[ 

2sinh(|^Bl) 



(0)|2 



(51) 



However, for our computation to be consistent, we must consider this expression in 
the large time limit. This is readily done in rescaling the areas A[ = A^/t and doing 
t — > oo. The final expression for the characteristic function of P({A^}) is quite simple: 

e i XUW) = (52) 



2 sinh 



B 'ii 



(B^ = £ J2i=i Bf ; when n — 1, (52) gives back Levy's result, eq.(5)). 

Owing to the form of B^, P({A'A-) will only be a function of the variable 



Yli=i(A' i } 2 (= A'), showing clearly that the (A^) 's are correlated. Its determination is 
reduced to the computation of the following integral [13]: 



/2n\ n / 2 1 r°° r n/2+i 

™-^')=(") J^rW^rfr (53) 

where </„ is a Bessel function. Closed form expressions can be given for odd n values. 
For instance: 

n = 3 P(A') = ^ tM *fr v^O (54) 

1 ; 2^cosh 2 (7r v / 3A') V ; 

r 5 tanh(7T V / 5A / ) - (71^^0(1 -3 tanh 2 (7r v / 5A / )) 

n = 5 P(A) = 4^ coshVv/5>) (55) 
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Now, if we consider the distribution of the sum of the areas A = J2?=i Ai, it is 
obtained by setting Bj = B, Vj. (52) leads to (A' = A/t): 

UBA<\ = (56) 

V I 2sinh(f) 1 ' 

With (5), we see that the sum of areas has, asymptotically, exactly the same distribution 
as the area enclosed by a single Brownian particle. In fact, we can compute [e %B - A ) for 
all t values (and not only when t — > oo). This is because, that time, the matrix B 
(By = B Sij) commutes with M. So, we are left with an { harmonic oscillator + 
uniform magnetic field } problem for each normal coordinate (except for Z\, that only 
feels a pure magnetic field). We get the result [14]: 



BA' 



B A Fi(B) 



2sinh(f)t L 2 Fi(0) 



A ^TnT ( 57 ) 



= 2^^rA~^^ (58) 




We recover (56) in the limit t -> oo (U?= 2 F i( B ) / F i(ty ~+ 1 when * -»■ oo). 



To close this section, it is interesting to consider the asymptotic law for the area A'j 
(= Aj/t) enclosed by a given monomer j. (52) gives: 

= r- ~ / B \ ( 60 ) 

' 20Isinh(^) 

It follows that Aj satisfies Levy's law (5) and scales like Remark that the area swept 
by the chain c.o.m., G, should scale like K On the other hand, for the same gaussian 
noise, we would get Aj ~ t if particle j was free (i.e. k — 0). The actual scaling of Aj 
is intermediate: this particle moves more freely than G but it is embedded in the chain, 
thus not completely free! 

Those considerations allow us to give a more precise sense to the statement: "the 
process is dominated by the c.o.m. motion". This one is true as long as we look at 
occupation times. However, when we study finer quantities like areas, this sentence 
must be corrected. Similar (even more dramatic) deviations will occur when we look at 
winding angles. 

To end up with areas, let us remark that the case of open trajectories can be treated 
exactly along the same lines as the one developed here, without additionnal difficulties 
(in particular, (48) still holds). We will not address this problem in the present work. 
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5 Winding angles distribution 



The last part of this paper will be devoted to the distribution P({9j}) (0j is the angle 
wound around O by particle j during a time t). We consider the same conditions as for 



Spitzer's law (z^°\ initial configuration, fixed, with Zj 7^ 0, Vj ; z, final configuration, 
unspecified ; t — > oo). 

We want to proceed as before and insert the constraint: 



lis [e s 

j=l \ 



2i Jo 



Z j Zj ^ i ^ ■ 



'3*3 



dr 



z j z j 



(61) 



in the Wiener measure (15). We are now faced with the problem of n harmonically bound 
particles submitted to the magnetic fields of n different point-like vortices located at the 
origin. The corresponding hamiltonian is: 



H Q + V X 



i=i 



v x = EmtS*-t 9 «)+E 



A? 



and the distribution P({9i}) is given by: 

dA 



P({0i}) 



TT iXjOj 

M 2 * 



i = l ^ZiZi 



dzdzF(z,z^) G x (z,z(°\t) 



(62) 
(63) 



G x (z,z<®,t) = {z 



(64) 



(65) 



Studying the limit t — > oo, we cannot develop directly as before a perturbation theory 
with V\\ this is because the last term in V x gives a divergent contribution [9]. Due to 



this term, all the eigenfunctions of H x must vanish in O at least as YY}=i \z. 
So, we redefine those eigenfunctions [9]: 



.IIAil 



The new hamiltonian acting on \l/ is 
Hx = H + V x 

VX(Z) : 



^ ( (A, - |A,|)-^ - (A, + |A,|)-^ 



i=l 



U(z)). 
(66) 

(67) 
(68) 



with a propagator G x 

G~ x (z,zM,t) 



-tH x 



U(z) 



G x (z,z^,t) 



(69) 
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(Go = Go). That time, the perturbation theory is properly defined and we can compute 
the characteristic function: 

G({A,}) ee (jZUw} = J dzdz (fl^oT^) F ^ z(0) ) Gl{z,z^\t) (70) 

with, symbolically, 

oo 

Gx= £(-l) m GS° (V x G^r (71) 

m=0 

Using integration by parts and also the relationship d Zi (jr^ = n8(zi), we first calcu- 
lated G({Aj}) up to 4 th order in V\, with the result: 

G({A,}) ~ e x ' 2 D(X) (72) 
D{X) = 1 + 

fn + l\(-X X 2 X 3 /3n 2 -l\ X 4 , n , oN \ 

+ hr hrr + -sr»-^r + ^r(3n 3 - 2n) - . . . (73) 



1! 2! 3! V 2 / 4! 

*=(l>*l) ln * ( 74 ) 

The prefactor e x l 2 comes out from U (z) in (70) when integrated over the final configu- 
ration: it will be present at all orders of the computation. Morover, (72) suggests that 
G({Aj}) is only a function of X: this is actually the case, as will be shown in the sequel. 

Let us consider the m th order term in (70,71) and suppose that we integrate, first, 
over z, z^ m \ z^™' 1 ^, . . . , z^ k+r) . Following the computation step by step, it is not difficult 
to convince oneself that the integration over z^ involves expressions of the form: 

J dz^dz^cf>(z^,T) V x (z^) G^ > (z < - k \z^ k ~ 1 \t k - t k -i) (75) 

z ( fe )|2 

where (p(z {k) ,T) = e — "~ e^E^^R T (76) 

and T = ti — t k , k + 1 < I < m. Let us call J k the result of (75). In the limit of long 
times, it reads: 



Jk — — \\ 



x 



\i=l 

2(4 — tk-i) 



x I — - oi-:^, t k - t fc _0 + - 1 0(z( fc - 1 ),T + 4 - t fc _i)V77) 

J + tfc — tfc_i / 



The m successive spatial integrations produce the factor (J2 |Ai|) m and, at the end, 
we are left with time integrals of the form: 

Em CKj f3^ 

I m (i m -i,...,i )(t)= dt m dt m -i... dti- r — — (78) 

J0 J0 JO [t;_, - t m _i . . . tii - h ti n 
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with Pi > and 

im-i = im-2 = . . . = i k = m ; i k _i = i k _ 2 = ■ ■ ■ = k = k ; = i t - 2 



We have proved, step by step, that; 



I m (im-i, ■ ■ ■ ,io){t) ~t- 



(Inty 



° nsHii-O 



ij = l; ... | 

(79) 



(80) 



Those considerations show that, actually, C({Aj}) is only a function of X. So, we 
can write D(X) = E™ =0 fl m I m , with a = 1 (see eq.(73)). 

Moreover, with the help of the above equation (80), and also looking at the tree 
structure exhibited in eq.(77), the following recursion relation can be shown: 



m—l 



y = 



n + 1 



It allows to write a closed form formula for D(X): 



-X/2 



D(X) = — 



y(e x - 1) cosh(X/2) + n sinh(X/2) 



(81) 
(82) 

(83) 



With (72) and, also, a rescaling of the angles = j^j) , we get the desired charac- 
teristic function: 



cosh(-u) + nsinh(-u) 



u 



i=i 



(84) 
(85) 



(with n — 1, we recover eq.(3)). 

We consider (84) as the main result of this paper. 

Finally, Fourier transformation shows that P({6'j}) is an "infinite sum of products 
of Spitzer's laws" (!) with highly correlated variables: 

I M 



fc=0 



'n - 1 
,n + 1 



n 



i 7T(2fc + 1) 

All the moments of this distribution are infinite (unless they trivially vanish). 



(86) 
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For a given particle j of the chain, we have: 



( 




1 



(87) 



cosh(Aj) + nsinh(|Aj|) 



that leads to: 




n(2k + 1) 



1 



1 + 




1 



) 



2 



(88) 



2k+l 



The difference with Spitzer's law is due to the presence of n in the denominator of (87). 

To shed some light on this problem, let us go back to the joint law (6) of small and 
big windings for the chain c.o.m.. What could we expect for the corresponding windings 
of particle j? With little effort, we can say that: 

i) The big windings will be roughly the same for both (when the chain is far from O, 
particle j follows the c.o.m. and winds around O in the same way). So, we keep 
A + unchanged in (6). 

ii) The small windings will be quite different. This is because particle j is artificially 
maintained in the vicinity of O: despite its higher mobility, it spends the same 
time as the c.o.m. in a given domain surrounding O. As a consequence, its small 
windings law will be broadened. Assuming that the remark following (7) holds, we 
get this broadening by changing |A_| into n|A_| in (6) (n is the ratio of the diffusion 
constants; of course, we don't say at all that (7) is the law of small windings!). 

Thus, our guess for particle j is: 



Setting A + = A_ = Xj, we recover (87). We are aware that this argument is strictly 
heuristic and that (89) remains to be proved. Nevertheless, we think that it allows to 
explain correctly the presence of n in (87). 

6 Conclusion 

Let us briefly summarize this work. We have computed explicitly the asymptotic joint 
laws of the occupation times, areas and winding angles of a chain of harmonically bound 
Brownian particles. 

For all these properties, we have shown that the scaling variables take the same 
general form as for the standard Brownian motion. However, a detailed study reveals 




(89) 
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important specific features that reflect a subtle interplay between the free c.o.m. motion 
- that strongly influences the whole chain properties - and the relative freedom of 
a given particle of the chain. For occupation times distributions, it appears that the 
c.o.m. satisfies the same law as a given monomer; now, for the areas, the scaling becomes 
slightly different and, finally, for the winding angles, the law itself is changed. Remark 
also that correlations are systematically present. 

Moreover, we observe that the scaling variables and the laws are very different from 
those met in our study of the attached Rouse chain [15] (in the latter case, 9 ~ t, A ~ y/i, 
and the winding angles are uncorrelated) . These differences are not so surprising since, 
in that case, we had no translation invariance. 



One of us (O.B.) acknowledges Dr. G. Oshanin for drawing his attention to this 
problem. 



16 



References 

[1] Rouse P E 1953 J. Chem. Phys. 21 1273 

[2] de Gennes P G 1979 Scaling Concepts in Polymer Physics (Ithaca) 

[3] Doi M and Edwards S F 1986 The Theory of Polymer Dynamics (Oxford) 

Grosberg A Y and Khokhlov A R 1994 Statistical Physics of Macromolecules (New- 
York, translated from russian) 

[4] Kallianpur G and Robbins H 1953 Proc. Nat. Acad. Sci. U.S.A. 39 525 

[5] Spitzer F 1958 Trans. Am. Math. Soc. 87 187 

[6] Levy P 1948 Processus Stochastiques et Mouvement Brownien (Paris) 

[7] Pitman J W and Yor M 1986 Ann. Prob. 14 733; 1989 17 965 

[8] Messulam P and Yor M 1982 J. Lond. Math. Soc. (2) 26 348 
Rudnick J and Hu Y 1987 J. Phys. A: Math. Gen. 20 4421 
Berger M A 1987 J. Phys. A: Math. Gen. 20 5949 
Brereton M G and Butler C 1987 J. Phys. A: Math. Gen. 20 3955 

[9] Comtet A, Desbois J and Monthus C 1993 J. Phys. A: Math. Gen. 26 5637 

[10] van Kampen N G 1981 Stochastic Processes in Physics and Chemistry (Amsterdam) 

[11] Riesz F and Sz-Nagy B 1955 Functional Analysis (Frederick Ungar Publ. Co.) 

[12] Feynman R P and Hibbs A R 1965 Quantum Mechanics and Path Integrals (New- 
York) 

[13] Hughes B D 1995 Random Walks and Random Environments (Oxford) 

[14] de Veigy A and Ouvry S 1992 Nucl. Phys. B 388 715 

[15] Benichou O and Desbois J 2000 submitted to J. Stat. Phys. 

e-mail: 

benichou@lptl.jussieu.fr 
desbois@ipno.in2p3.fr 



17 



